Are gravitational waves from giant magnetar flares observable? 
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Are giant flares in magnetars viable sources of gravitational radiation? Few theoretical studies 
have been concerned with this problem, with the small number using either highly idealized models 
or assuming a magnetic field orders of magnitude beyond what is supported by observations. We 
perform nonlinear general-relativistic magnetohydrodynamics simulations of large-scale hydromag- 
netic instabilities in magnetar models. We utilise these models to find gravitational wave emissions 
over a wide range of energies, from 10^" to 10*'' erg. This allows us to derive a systematic relation- 
ship between the surface field strength and the gravitational wave strain, which we find to be highly 
nonlinear. In particular, for typical magnetar fields of a few times 10^^ G, we conclude that a direct 
observation of /-modes excited by global magnetic field reconfigurations is unlikely with present or 
near-future gravitational wave observatories, though we also discuss the possibility that modes in a 
low-frequency band up to 100 Hz could be sufficiently excited to be relevant for observation. 
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Introduction. Bursts and occasionally giant flares in 
soft-gamma repeaters (SGRs) and anomalous X-ray pul- 
sars (AXPs) are commonly understood as consequences 
of magnetic field activity in magnetars, neutron stars en- 
dowed with strong magnetic fields Jll. These violent 
events have been considered in the literature as pos- 
sible efficient sources of gravitational radiation, since 
they are very compact, and since the high electromag- 
netic luminosity may represent only a small part of the 
overall energy contained in the mechanism. It is thus 
not surprising that operational gravitational wave ob- 
servatories, in particular LIGO and VIRGO, have been 
employed to establish upper limits on SGR bursts and 
storms. Presently, the best empirically derived upper 
limits are 1.4 x 10^^ erg for the /-mode, and 3.5 x 10^^ erg 
for a white noise band around 100 — 200 Hz 121 • For 
the /-mode, gravitational wave emission of this mag- 
nitude would imply a very substantial excitation of the 
star. Can luminosities of this order actually be realized 
in a nearby giant flare? 

There has been little theoretical work addressing this 
question, for a large part because the system is very 
complex and the giant flare mechanism is not yet well 
understood. Possible flare mechanisms are discussed in 
||3HZl and fall into two broad categories: internal mecha- 
nisms based on a large-scale rearrangement of the inte- 
rior magnetic field, and external mechanisms likely fol- 
lowing magnetic reconnection in the crust and magneto- 
sphere. In this work, we will assume a scenario involv- 
ing a global restructuring of the magnetic field. 

Two direct pieces of evidence correlate the electro- 
magnetic signal from magnetar flares with their mag- 
netic field dynamics - energy budgets associated with 
the flare luminosities (S] and observations of quasi- 
periodic oscillations in the tails of the flares IISHTOl. 
An optimistic estimate for gravitational wave emissions 
from magnetic field dynamics can therefore be attained 
by assuming a large-scale, dynamical rearrangement of 



the core magnetic field inside the star. loka IfTTJI has in- 
vestigated the maximum gravitational wave energy re- 
leased by the change in moment of inertia induced by 
such a mechanism and placed an upper limit of about 
10*^ erg under ideal conditions, including optimistic val- 
ues of the internal magnetic field. More recently, Corsi 
and Owen [12 1 found similar values to be possible under 
more generic conditions, still tapping into the full en- 
ergy reservoir associated with an instantaneous change 
in the magnetic potential energy of the star. In contrast. 
Levin and van Hoven [7J do not find /-mode detection 
to be very likely in the near future. Their model is based 
on the aforementioned external mechanism triggering 
/-modes in the star. In the present paper, we focus on 
the other side of the problem - attempting to trigger 
large-scale mass motions that could generate gravita- 
tional radiation through a global rearrangement of the 
internal magnetic field. 

Our recent publication on dynamical instabilities ||T3| 
has focused on the dynamics and evolution of mag- 
netic fields inside relativistic stellar models, and we 
have followed the development of these instabilities un- 
til saturation and ring-down to establish the nature of 
the quasi-stationary states which result from these pro- 
cesses. Most recently, Ciolfi et al [14| have performed 
numerical simulations of the same hydromagnetic insta- 
bilities, and concluded that giant flares could give rise 
to observable gravitational radiation, employing a stel- 
lar model with a surface field strength of 6.5 x 10^^ G 
to reach this conclusion. Actual SGR field strengths are 
closer to 1/60 of this value |1 15 1, and the luminosity 
depends on the field strength in a highly nonlinear fash- 
ion, as we shall see below. 

In this paper, we will address the following ques- 
tions: (i) Assuming a large-scale rearrangement of the 
core magnetic field to be involved in a giant flare, how 
does the gravitational wave luminosity scale with the 
field strengths? (ii) As a consequence of this relation. 



can we expect to observe magnetar giant flares in gravi- 
tational wave detectors? 

Model. In order to investigate gravitational wave lu- 
minosities, we consider a large-scale restructuring of an 
initially purely poloidal magnetic field. While an ac- 
tual magnetar should be quasi-stationary before the flar- 
ing event, we employ a hydromagnetically unstable star 
to investigate the gravitational wave signal associated 
with the instability. This provides us with a simplified 
model that mimics a catastrophic reconfiguring of the 
internal magnetic field during a magnetar flare. 

As with our recent publication [13J, we are using 
the graphics processing unit (GPU)-accelerated HORI- 
ZON code iri6l for general relativistic magnetohydrody- 
namics in the Cowling approximation, in conjunction 
with the LORENE library for the construction of mag- 
netized equilibrium neutron stars IITtI . HORIZON is 
based on the Thor code flF. "191, but employs CPUs 
for high-throughput parallel processing. Gravitational 
waves are extracted using the first-moment form of the 
quadrupole formula. For more details on the numeri- 
cal method and stellar models see Lasky et al. |13|, Zrnk 

Our initial models are relativistic equilibrium n = 1 
polytropes (mass 1.3 M©, radius 15 km) with purely 
poloidal magnetic fields of varying strength. In partic- 
ular, we have investigated a sequence of models with 
polar surface magnetic field amplitudes between B„gif, = 

3.1 X 10^^ G and 5.5 x 10^^ G corresponding to magnetic 
fleld strengths of 1.6 x 10^^ G to 2.7 x 10^^ G in the center 
of the star While typical magnetar field strengths reside 
near the low end of this spectrum, we have decided on 
a sequence of models to gain fundamental insights into 
the relation of initial field configuration and the gravita- 
tional wave amplitude. 

Purely poloidal fields are known to be d5n-iamically 
unstable 11131 I20L and we follow the evolution through- 
out development, saturation and ring-down of the insta- 
bility This is particularly challenging for models with 
lower field strengths, since we need to follow the sys- 
tem for many Alfven times in all cases. 

Results. All our models exhibit the poloidal field in- 
stability as expected. In figure 111 we show the measured 
gravitational wave strain /ix at lOkpc in a particular 
model with B„gig = 1.8 x 10^^ G. Other components and 
other models in the sequence show a similar structure, 
although the growth timescale of the instability and the 
strain amplitude are different. An analysis of the sig- 
nal spectrum shows that the fundamental quadrupole 
/-mode is excited by the magnetic field instability. In- 
triguingly we also see some evidence that low-frequency 
modes in the range of 100 Hz may contribute to the grav- 
itational wave signal. This evidence is not yet conclusive 
and warrants further investigation, in particular very 
long-term simulations in the order of a second or more 
to gain a sufficiently high resolution in this part of the 
spectrum. 

Figurel2]collects the gravitational wave amplitudes as 
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Figure 1: Gravitational wave signal obtained from one of our 
simulations. This particular stellar model has an initial sur- 
face magnetic field strength of 1.8 x 10^^ G and develops a 
poloidal field instability during the first 50 ms (see also |13|) 
which restructures the global magnetic field inside the star. 
The gravitational wave signal, here h x for an assumed source 
distance of lOkpc, has the expected 1.8 kHz oscillations asso- 
ciated with the /-mode, and also exhibits low frequency com- 
ponents which are further discussed in the text. 



a function of the magnetic field strength. The approx- 
imate values and error bars for the strains are derived 
from the time variation in the h x signals (w^hich are simi- 
lar in amplitude to h^) after the instability has saturated. 
We find an approximate power-law relation between h 
and Bpoig given by 
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Most of the energy in the signal is in the /-mode. As- 
suming a gravitational wave damping time of approx- 
imately 100 ms (e.g. ||2T1|22| ), we find a corresponding 
power-law relation for the energy emitted in gravita- 
tional radiation: 
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From this result, which can be directly compared with 
observations [2], we make three notes: (i) The ratio of 
gravitational wave energy extracted from our simula- 
tions to the electromagnetic energy emitted in a flare fT) 
is low, in the order of Egm/Eem ~ lO''^ erg/ 10** erg = 
10~^. (ii) The gravitational wave amplitude is a 
highly nonlinear function of the surface magnetic field 
strength, (iii) Typical magnetar field strengths of 10 G 
give rise to strains below lO^'^^ and energies lower than 
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Figure 2: Gravitational wave strain in relation to surface mag- 
netic field strength at the stellar pole. This plot represent re- 
sults from our entire sequence of simulations, with surface 
magnetic fields ranging from 3.1 x 10^^ G to 5.5 x 10^^ G. The 
signal amplitude is found to depend strongly on the magnetic 
field strength, ranging from almost 10^^^ for our most extreme 
model down to less than 10^^^ for models closer to realistic 
magnetar strengths. The dashed line is a power-law fit to the 
data (see text). 



10**^ erg for a source at 10 kpc, even if we assume a catas- 
trophic global restructurmg of the field to be associated 
with a giant flare. 

We have indicated the signal-to-noise ratios we derive 
from our data for different detector sensitivity curves 
in figure [3] Here we have plotted the signal ampli- 
tude VT \h{f)\, where T is the damping time of the 
oscillation and h{f) is the Fourier transform of hy, as 
a function of the frequency for various magnetic field 
strength models. At approximately 1.8 kHz we show the 
signal amplitude for the /-mode as excited by the hy- 
dromagnetic instability and assuming a damping time 
of 50 ms < T < 200 ms. In the lower part of the 
spectrum we plot the other maximal mode seen in the 
Fourier transform, assuming a damping time between 
of 10 ms < 7" < 1 s. Over this we plot the entire 
spectrum (assuming T = 100 ms) for two models with 
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tively. Finally, we also plot the root of the noise power 
spectral density, \/|S;,(/)|, as a function of frequency for 
the LIGO, AdvLlGO and ET detectors Il23l . 

The amplitude signal-to-noise ratio, defined by 

a/T \h{f)\ I I Sft (/) I can be read off the graph as the 
ratio between the signal and the noise curve for the re- 
spective detector. The main conclusion from this figure 
is this: Assuming that a giant flare is associated with a catas- 
trophic large-scale rearrangement of the core magnetic field, 
the gravitational wave signal associated with f -modes are not 
observable with present or near-future gravitational wave ob- 
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Figure 3: Signal amplitude v'T|I!(/)| against the oscillation 
frequency. The colored boxes on the right represent the maxi- 
mum for the /-mode assuming a constant periodic source last- 
ing between 50 and 200 ms. The colored boxes on the left rep- 
resent the maximum mode seen in the Fourier transform for 
any given frequency below the /-mode frequency assuming a 
constant periodic source lasting between 10 ms and 1 s. These 
scale with \/T, implying one can easily extrapolate to alter- 
native values of the damping time. We have further plotted 
the entire spectrum (assuming a damping time of T = lOO ms) 
for both the Bp,,;,, = 1.8 x 10^'' G model (blue line) and the 

^pole = 8-8 ^ 10^^ G model (black line). The noise power spec- 
tral density curves have been taken from the review article of 
Sathyaprakash and Schutz Ii23j . 



servatories. 

This statement is in line with the conclusion of |7| 
concerning the efficiency of the internal mechanism in 
exciting /-mode radiation. Moreover, it is in contrast 
with the scenarios presented in |11 , 12J, which are both 
based on equilibrium models as opposed to full dynam- 
ical simulations. While we cannot exclude that differ- 
ent initial model configurations could produce higher 
levels of /-mode radiation, we consider it, at present, 
unlikely that the result will change by many orders of 
magnitude. 

Our results can also be compared with the numeri- 
cal study of [14J which obtains a much more optimistic 
prediction than we do, although they investigate the 
same kind of instability. However, the magnetar model 
considered by these authors has a very strong surface 
magnetic field of 6.5 x 10^^ G, which is even stronger 
than the most extreme case we have considered here. 
Such a high magnetic field strength favors a coupling of 
Alfven modes into /-modes because the Alfven speed 
approaches the speed of sound. Clearly, the observed 
soft-gamma repeaters do not fall into this category. 

Discussion. We have performed three-dimensional 
general-relativistic magnetohydrodjn-iamics simulations 
of neutron stars, inducing large-scale reconfigurations 
of the internal magnetic field to model gravitational 
wave emissions from magnetar flares. We have found 



power-law relations governing the surface, polar mag- 
netic field strength as a function of the maximal grav- 
itational wave strain (equation [T| and also the energy 
emitted in gravitational waves (equation l2l. We find 
that the gravitational wave emissions due to /-mode ex- 
citations are unlikely to be observed in current or near- 
future gravitational wave observatories (figure l3|. 

There are two factors which could modify our above 
conclusion. The first is that other equations of state 
and initial magnetic field topologies could be investi- 
gated, e.g. those with a dominating toroidal field com- 
ponent. Mixtures of poloidal and toroidal fields in neu- 
tron stars may be stable for a range of energy ratios 
Il24l and could indeed represent typical field equilibria. 
However, there is a possibility that configurations with a 
strong toroidal component (and comparatively weak sur- 
face dipole field) could be unstable and lead to higher 
coupling into the /-modes and consequently stronger 
gravitational wave luminosities, w^hile still being consis- 
tent with SGR spindown rates. Recent studies of hydro- 
magnetic equilibria in stratified neutron stars [25) i26J , 
however, cast a shadow of doubt over this possibility. 

The other factor is the possible low-frequency compo- 
nent of the signal (see also [,7, .27J ). As discussed above. 



we have found indications (but no solid confirmation) 
of spectral components in the band typically associated 
with g-modes or Alfven modes, i.e. below 200 Hz. This 
is particularly interesting because the detectors are most 
sensitive in this regime. Moreover, these modes have 
much longer damping times and are assumed to be part 
of the post-flare signal. We leave these possibilities for 
future work. 
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